Gravitational radiation and angular momentum flux from a slow rotating dynamical black hole 
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A four-dimensional asymptotic expansion scheme is used to study the next order effects of the nonlinearity 
near a spinning dynamical black hole. The angular momentum flux and energy flux formula are then obtained 
by constructing the reference frame in terms of the compatible constant spinors and the compatibility of the 
coupling leading order Newman-Penrose equations. By using the slow rotation and small-tide approximation 
for a spinning black hole, we chose the horizon cross-section is spherical symmetric. It turns out the flux formula 
is rather simple and can be compared with the known results. Directly from the energy flux formula of the slow 
rotating dynamical horizon, we find that the physically reasonable condition on the positivity of the gravitational 
energy flux yields that the shear will monotonically decrease with time. Thus a slow rotating dynamical horizon 
will asymptotically approaches an isolated horizon during late time. 
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I. INTRODUCTION 



Null infinity and black hole horizon have similar geometrical 
properties. They are both three dimensional hypersurfaces and 
have the gravitational flux across them. The physical proper- 
ties of null infinity can be studied in the conformal spacetime 
with finite boundary. Thus the conformal method provides an 
alternative way to study Bondi-Sachs gravitational radiation 
near null infinity, which was first proposed by Penrose lfl9ll . 
The boundary of a black hole is asymptotically non-flat and 
one may not be able to apply the conformal method to study 
the boundary problem of a dynamical black hole. Rather than 
using the symmetry for the whole space-time to locate the 
boundary of a black hole, Ashtekar et al. use a rather mild 
condition on the symmetry of the three dimensional horizon 
01 Si- This quasi-local definition for the black hole boundary 
makes it possible to study the gravitational radiation and the 
time evolution of the black hole. 

In this paper, we use the Bondi-type coordinates to write 
the null tetrad for a spinning dynamical horizon (DH). The 
boundary conditions for the quasi-local horizons can be ex- 
pressed in terms of Newman-Penrose (NP) coefficients from 
the Ashtekar's definition on DH. Unlike Ashtekar et al's (HIH 
three dimensional analysis, we adopts a 4-dimensional asymp- 
totic expansion to study the neighborhoods of generic iso- 
lated horizons (IHs) and dynamical horizons (DHs). Since 
the asymptotic expansion has been used to study gravitational 
radiations near the null infinity H O], it offers a useful 
scheme to analyze gravitational radiations approaching an- 
other boundary of space-time, black hole horizons. We first 
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set up a null frame with the proper gauge choices near quasi- 
local horizons and then expand Newman-Penrose (NP) coeffi- 
cients, Weyl, and Ricci curvature with respect to radius. Their 
fall-off can be determined from NP equations, Bianchi equa- 
tions, and exact solutions, e.g., the Vaidya solution. This ap- 
proach allows one to see the next order contributions from the 
nonlinearity of the full theory for the quasi-local horizons. 

We have shown that the quasi-local energy-momentum flux 
formula for a non-rotating DH by using asymptotic expansion 
yields the same result as Ashtekar- Krishman flux IU41IT5I1 . For 
slow rotating DH, we have presented our results in 11411 . how- 
ever, we use an assumption of vanishing NP coefficient A on 
DH. Furthermore, the energy-momentum flux formula has a 
shear (NP coefficient er) and a angular momentum (NP coef- 
ficient 7r) coupling term. Since it is unclear whether the exis- 
tence of this term carries any physical meaning or it may due 
to our assumptions, we thereby extend our previous work on 
IHs and DHs into a more general case. 

An algebraically general structure (Petrov type I) of spacetime 
is thought to be related with gravitational radiation for an iso- 
lated source and can tell us more about the inner structure of 
the gravitating source. The Weyl scalars ^k,k = 0, .., 4 can 
be expanded in terms of an affine parameter r along each out- 
going null geodesic based on assumption of compatification 
of null infinity fKH]. Here = 0(r fe ~ 5 ),fc = 0, ..,4 and 
one may find that it peels off more and more when moving 
inward along null ray. From Ashtekar's definition of isolated 
horizon, it implies that v&Oj^i vanishes on horizon. There- 
fore the space-time is algebraically special on isolated hori- 
zon. However, spacetime may not be algebraically special for 
an arbitrary DH. The corresponding peeling theorem for an 
arbitrary DH is crucial for our gravitational radiation study. 
Due to the difficulties of knowing the fall-off of Weyl scalars, 
we use Kerr- Vaidya solution to serve as our basis for choosing 
the fall-off of Weyl scalars, which is ^q, vanishing on DH, 
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in our previous work on slow -rotating DH [14']- So space- 
time structure on slow-rotating DH is still assumed to be al- 
gebraically special. However, according to the gravitational 
plane-wave solutions, ^o an d ^4 indicate the ingoing and out- 
going gravitational waves, respectively. It seems physically 
unsatisfactory to assume $o vanishing on DH. Moreover, the 
algebraically general space-time allows four roots of the equa- 
tion, which correspond to the principal null-directions of Weyl 
scalars, and describes the gravitational radiation near the grav- 
itating source. Therefore it would be more reasonable for one 
to consider an algebraically general space-time on an evolving 
DH. From the reduction and the decoupling of the equations 
governing the Weyl scalars, instead of assuming ^i = 
on DH, we set "Ji, ^3 vanishing on a spinning DH. This is a 
similar setting with the perturbation method (Also see Chan- 
drasekhar |01). 

We present the results of asymptotic expansion for a spinning 
DH in Sec. [HI] However, it maybe too general to yield some 
interesting physical results. By considering the small-tide and 
slow rotate of DH and using slow rotate Kerr solution as a ba- 
sis, we use two sphere conditions of DH cross section for our 
later calculation. The NP coefficient Ao (shear for the incom- 
ing null tetrad n) on DH is no longer assumed to be vanished 
when calculating the flux formula. The index on NP coef- 
ficients denotes their values on DH. Directly from non-radial 
NP equations, we find that ao and tto coupling terms can be 
transformed into ttq terms only, so the problems of our previ- 
ous work fllll are resolved. 

Though the exact solution for a stationary rotating black hole 
has been found near fifty years, the spacetime with rotation 
remains its ambiguity and difficulty for quasi-local mass ex- 
pressions and boundary condition. For example, the existence 
of angular momentum will not change the boundary condi- 
tion for the null infinity, however, it will affect the boundary 
condition of a black hole. Among the well-known quasi-local 
mass expressions named Komar, Brown- York and Dougan- 
Mason, only Komar integral of the quasi-local mass for an 
arbitrary closed two surface can go back to the unique New- 
tonian quasi-local mass lfl6ll . Unfortunately, Brown- York and 
Dougan-Mason mass can return to the unique surface integra- 
tion of the Newtonian mass in the covariant Newtonian space- 
time only for the spherically symmetry sources. In GR, quasi- 
local mass expressions for Kerr solution disagree one another 
JH]. Different quasi-local expressions give different values of 
quasi-local mass for Kerr black hole. At null infinity, there is 
no generally accepted definition for angular momentum lfl7ll . 
Unfortunately, no explicit expression for Bramson's angular 
momentum in terms of the Kerr parameters m and a is given 
@]. We use Komar integral to calculate angular momentum 
since it gives exactly ma for Kerr solution. Although differ- 
ent quasi-local expressions yield the different results for Kerr 
solution, our main motivation is to analyze and discuss the 
compatibility of the coupling NP equations from asymptotic 
expansions. We both calculate quasi-local mass and flux for a 
spinning DH based on two spinors (Dougan-Mason) and Ko- 
mar integral. It is found that these two expressions yield the 
same result. 



Bondi and Sachs use no-incoming radiation condition for 
gravitational wave on null infinity iH 0]. However, no- 
incoming radiation condition is only true for linearized the- 
ories, e.g., electrodynamics and linearized GR, as to exclude 
the incoming rays. The incoming pulse waves do not destroy 
the asymptotic conditions for null infinity since they are ad- 
mitted by formalism. Their existence may play an important 
role in the interpretation of the new conserved quantities (NP 
constants) IIllll2tl . The interpretation and physical meaning 
of these constants have been a source of debate and contro- 
versy until today. Some physical discussions and application 
of them can be found in |@, QjJ. Despite the vagueness of 
the physical meaning of these conserved quantities, in the full 
nonlinear gravitational theory, the mass and momentum are no 
longer absolutely conserved and can be carried away by the 
outgoing gravitational wave, so as to give a positive energy 
flux at infinity. Here we consider a spacetime inner bound- 
ary, e.g, a spinning DH in this paper. By the aid of using 
asymptotic constant spinor to define spin frame as the refer- 
ence frame for our observation, mass and angular momentum 
flux can be calculated. According to the coupling NP equa- 
tions from the asymptotic expansion analysis, such a system 
will gain energy and will cause the radius of the black hole 
to increase. From similar argument, the outgoing waves do 
not change the boundary conditions of the quasi-local hori- 
zons (DHs or IHs) and make no contribution to flux, while 
incoming wave will cross into DHs. The existence of incom- 
ing wave indicates the difference between IHs and DHs. The 
mass and momentum are carried into black hole by the incom- 
ing gravitational wave. 

In the Ashtekar-Kirshnan's 3-dimensional analysis, it gives no 
constraints on co's time evolution. However, through our 4- 
dimensional asymptotical expansion scheme, we observe that 
the Co will monotonic decrease with time once the positivity 
of gravitational energy flux is hold on a slow rotating DH. It 
means that the slow rotating DHs will gradually settle down 
to IHs as Co approach zero, which is physically reasonable. It 
is similar to a physical assumption saying that the mass loss 
cannot be infinite large for null infinity 11311 . However, rather 
than assuming that mass loss cannot be infinite large we ob- 
tain this result directly from asymptotic expansion analysis for 
a slow rotating DH together with the physical arguments. Fur- 
ther from the commutation relations, we find that the horizon 
radius of a slow rotating DH will not accelerate. The slow 
rotating dynamical horizon increases with a constant speed. 
There is one more interesting point about the peeling proper- 
ties for a slow rotating DH. It is known that the peeling proper- 
ties refer to different physical asymptotic boundary conditions 
of a slow rotating black hole. By comparing our current work 
to previous one IU5I1 . which have different peeling properties, 
and also due to the monotonic decrease of ao, we propose that 
the setting of Weyl scalars in this work makes it excludes the 
possibility to absorb the gravitational radiation from near by 
gravitating source. 

The plan of this paper is as follows. In Sec. [TTJ we re- 
view the definition of DH and express Ashtekar-Krishnan's 
3-dimensional analysis of DH in terms of NP coefficients. The 
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gauge choices and boundary conditions of a spinning DH are 
applied to the asymptotic expansion in Sec. [TTT] In Sec. |IV| 
we first examine the gauge conditions of slow rotating Kerr 
solution in Subsec. IIV Al Later we use the two sphere condi- 
tion for a slow rotating DH with small tide in Subsec. IIV Bl 
The results of asymptotic expansion are largely simplified by 
considering DH's cross section as a two sphere. Angular mo- 
mentum and its flux for a slow rotating DH are calculated by 
using Komar integral in Sec. [V] Energy-momentum and its 
flux of a slow rotating DH are obtained in Sec. [VTJ We first 
calculate mass and mass flux by using Komar integral in Sub- 
sec. I VI Al Then, mass and mass flux of a slow rotating DH is 
calculated by using two spinor method in Subsec. IVIBI The 
time evolution of shear flux and its monotonic decrease is dis- 
cussed here. We find that either Komar integral or two spinor 
method yields the same result. 

In this paper, we adopt the same notation as in U 0] for de- 
scribing generic IHs and DHs. However, we choose the differ- 
ent convention (H ), which is a standard convention for 

the NP formalism lfl8ll . The necessary equations, i.e., com- 
mutation relations, NP equations and Bianchi identities, for 
asymptotical expansion analysis can be found in p. 45-p. 51 
of 0]. We use "=" to represent quantities on a dynamical 
horizon (ignore 0(r')) and use "=" to represent quantities on 
a slow rotating horizon (ignore 0(a 2 )). 



II. ASHTEKAR DYNAMICAL HORIZON 

A. The dynamical horizon 

The generic IHs are taken as the equilibrium state of the DHs. 
The DHs can be foliated by marginally trapped surface S. 
Therefore, the expansion of the outgoing tetrad vanishes on 
DHs. 



a. Definition A smooth, 3 dimensional, space-like sub- 
manifold H of space-time is said to be a dynamical horizon 
(DH) if it can be foliated by a family of closed 2-manifold 
such that: (1) on each leaf, S, the expansion of one null 
normal £ a vanishes, (2) the expansion &( n ) of the other null 
normal n a is negative. 



From this definition, it basically tells us that a dynamical hori- 
zon is a space-like hypersurface, which is foliated by closed, 
marginally trapped two surface. The requirement of the ex- 
pansion of the incoming null normal is strictly negative since 
we want to study a black hole (future horizon) rather than a 
white hole. Also, it implies 

Re[p]=0,ReM < 0. (1) 



B. Dynamical horizon in terms of Newman-Penrose 
coefficients 

If we contract the stress energy tensor with a time-like vector, 
then in components T k o represents the energy flux of matter 
field. Therefore we can use a time-like vector T a and contract 
it with the stress energy tensor to define the flux of the mat- 
ter energy. Here we are more interested in the energy of the 
matter field associated with a null direction. One can thus cal- 
culate the flux of energy associated with £ a = N£ a . The flux 
of matter energy across H along the direction of I is given by 

^matter \= [ T ab T a ^d 3 V. (2) 
J H 

The dynamical horizon is a space-like surface, the Cauchy 
data ( ( 3 'q a b, K ab ) on the dynamical horizon must satisfy the 
scalar and vector constraints 

H s : = ^R + K 2 -K ah K ab = 16TrT ab T a T\ (3) 
: = V b (K ab - K ( 3 Y 6 ) 
= V b P ab = 8irT bc T c (3) <f 6 (4) 

where P ab := K ab - K ^q ab . 

If the dominate energy condition is satisfied, it turns out that 
H has to be a space-like hypersurface [2J|. The unit time- 
like vector that is normal to H is denoted by T a and the unit 
space-like vector that orthogonal to the two sphere and tangent 
to H is denoted by R a . In order to study them in terms of 
Newman-Penrose quantities, they can be defined by using the 
null normals l a and n a . Therefore, 

T a = -^=(£ a + n a ), R a = -=(l a - n a ) (5) 

where T a T a — 1, R a R a = —1. The four metric has the form 

g ab = n a e b +e a n b + (2) q ab (g) 

= T a T b -R a R b + (2 Y b . (7) 

The three metric ^q ab that is intrinsic to dynamical horizon 
His 

(3) q ab = g ab _ T a T b = (2) q ab _ R a R b (g) 

The two metric Wq ab that is intrinsic to the cross section two 
sphere S is 

(2) q ab = (3) g ab + R a R b = _( m ^b + ^ m &). (Q) 

The induced covariant derivative on H can be defined in terms 
of 4-dimensional covariant derivative V a by 

V b V a := {3) q b C (3) <Za d V c F d , (10) 

so the three dimensional Ricci identity is then given by 

(3) iUc d w d = -[V a ,V h ]w c . (11) 
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The induced covariant derivative on cross section S can be 
defined in terms of 4-dimensional covariant derivative V a by 

®V b V a := ( 2 V (2 V d V c F d . (12) 
The extrinsic three curvature K ab on H is 

K ab = {3) q(a c ( % d V c T d = V a T b - T a ^a b 

where < 3 )a h = T c V c T h . One can also introduce the extrinsic 
two curvature ^K ab on S by 

{2) K ab = Wq (a c Wq b ) d V c R d = V a R b + R a ^a b , 

where (2) a 6 = R c V c R b = R (c ^q b d) V c R d . After a straight- 
forward but tedious calculation, we can write the extrinsic cur- 
vature in terms of NP spin coefficients. Here we present the 
general extrinsic three curvature and two curvature without 
any assumption of gauge conditions. The extrinsic three cur- 
vature K ab is 

K ab = (3 V c V c T b) = Wq {a c V c T b) -R {a R c V c T b) 
= A^qab + S ab + 2W {a R b) +BR a R b (13) 

where 



S a b = — 7=[(cr - X)m a m b + C.C.] 
V2 



.4 



Rep — Re/z 

W a := - {2 ka C KcbR b 

= v - t - 7T - 2(a + ~j3)]m a + C.C.. 



B = -V2(Ree-Re7), 

and C.C. denotes the complex conjugate terms. The extrinsic 
two curvature ( 2 'K ab is 

(2) K ab = Wq {a c V c R b) = W q (a c{3 hc d(3) q b ) e V a R e 

= {2 ha d{2) q b e V d Re = \ {2) K^qab + (2) S ab 



where 



(2 



(14) 



{2) S ab = -=(W + X)m a mb + C.C. (15) 
v2 

The calculation of two acceleration ( 2 'a a yields 

Wa a = R b V b R a = R^ c ^q a d ^V c R d = Cm a + Cm a , 
where 



C = - v + n - r), 

so the two acceleration is tangent to S. 



(16) 



We now perform 2+1 decomposition to study the various 
quantities on H. The curvature tensor intrinsic to S is given 
by 

-^R abc d = ~^q a f^q b ^ 2 \ c k ^q/^R fg ^ 
-^K ac ^K b d + ^K bc ^K a d , 

which is the Gauss-Codacci equation. This leads to the rela- 
tion between the scalar three curvature ^R and the scalar two 
curvature ( 2 >R 

_(3) R = _(2) R _(2) R 2 + (2) K V) K ab _ 2 V a a a (17) 

where a a := R b V h R a - R a V b R b = ( 2 V - R a ^K. 
From (1171 , we obtain the Einstein tensor on H 

- 2 ^G ab R a R b = -^R +< 2 > K 2 - W K ab ^K ab . (18) 

The expansion of the out going tetrad £ a can be calculated to 
yield 

Q w :=-l( Jt + ji) = J-[K+ &K + B). (19) 

Now we use the following relations (f20b-(f23T> to calculate 
Hg + 2R a H v , where H$ and Hy are the scalar and vector 
constraints defined in (0) and 

K = 2A — B, (20) 
{2) K = -K - B + 2\/29(*) = -2A + 2-^/29^,(21) 
K ab K ab = 2A 2 + S ab S ab -2W a W a + B 2 , (22) 

( 2 ^ a() < 2 )/v ab = - ^K 2 + ^S ab ^S ab . (23) 

From the momentum constraint and use integration by 
parts, we get 

R b V a P ab = V a p a - P ab V a R b (24) 

where /3 a := K ab R b - KR a . Thus, 

7 a := a a + j3 a = R b V b R a - W a - 2V29 (0 i?°. (25) 

For a general space-time, the matter energy flux can be calcu- 
lated as following 

H s + 2R a H v = {Z) R + K 2 - K ab K ab + 2R a V b P ab 
= MR +< 2 > K 2 -< 2 > Kab ^K ab + K 2 - K ab K ab 

-2P ab V a R b + 2V al a (Use 03) 
= ^ R - a ab a ab + 2W a W a -2W a(2) a a 

+AQ {l) {Q {l) ~V2B) + 2V al a - 

By applying the 2+1 decomposition on the covariant deriva- 
tive T> and using integration by parts, we have 

2V al a = 2V a ( (2 V - W a - 2^/2Q (l) R a ) 
= 2( (2 V ^a a - W a < 2 >a a - V a Q (l) R a + \ ( 2) ^6 W ), 
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where the term ^V a (^ a a — W a ) has been discarded since it 
will vanish due to the integration over a compact two surface 
S, and then 

2(W a W a - W a ^a a +V al a ) 
= 2(W a - ^a a )(W a - Wa a )~2V a e (e) R a - (2) KQ (e) . 

Here we can define 

C a := W a - Wa a = -V2W q i d R c) V c £ d 

= \[K-t- {a + P)]m a + C.C. (26) 

Finally, we get 

H s + 2R a H v = «R- <r ab <r ab + 2( a ( a 

- 2V a Q {l) R a + 9 W (- {2) K + 46 w - 4V2B) 

where 

(Tab = A=(S a b + {2) Sab) = V2am a m b + C.C. 
v2 

is the shear of null normal l a . This equation is completely 
general. On the dynamical horizon, the outgoing expansion 
Q(f, vanishes. It then becomes 1 

F matter =-i- / N( a ab a ab + 2( a C)d 3 V. (27) 

Iott Jam 

If the gauge condition 

K=7T-T=7T- (a + /3)=0 (28) 

is satisfying, where = denotes the equating on DH, then £ a = 
—Trm a + C.C. and £ a £ a = — 2irW. So the flux formula in 
terms of NP in this gauge is 
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Fmatte r =T^- / N( {2) R - Aaa - 47T7f) d 3 V. (29) 

lOTT J AH 



C. Angular momentum flux and energy fluxes 



By contracting the vector constraint H v with the rotational 
vector field ip a , which is tangential to S, we can obtain angu- 
lar momentum of a black hole. Then we integrate the result- 
ing equation over the region of AH and use the integration by 
parts together with the identity C$ ^q a b = 2V^ a ijj b ). It leads 
to 



dJ — Js 1 — Js 2 



— I K ab ip a R b dS — -~- (j) K ab ip a R b dS 



Sir 



S-2 



8- 



Si 



/ (T ab T a ^ b + T^-P a % {3) dab)d 3 V. (30) 
J ah !6tt 



1 Our expression has some minus sign different from Ashtekar's expression 
because of convention. 



The angular momentum associated with cross-section S is 

= — §tt $s Kabit> a R h dS where ip a need not be an axial 
Killing field. The flux of angular momentum due to matter 



fields -Fitter and gravitational waves Fg TSW are 



x m 



T ab T a ^ b d 3 V, 



(31) 



AH 



= -TT I {3) 1abd 3 V, (32) 

■LO 71 " J AH 

and the balance equation Jg — Jlj i — Flatter + Fgav, which 
describes the difference of angular momentum between two 
cross section, is due to the matter radiation and gravitational 
radiation. 

Each time evolution vector t a defines a horizon energy E A . 
From equation d27l i. we find the total energy flux is the com- 
bination of the matter flux and gravitational flux 

^matter + ^grav = T~ f N {2) R d?V (33) 
16tt J ah 

where the matter flux is equation and the gravitational flux 
is 

F e rav = ^ / N(a ab a ab - 2( a C)d 3 V. (34) 
16tt J ah 



(35) 



If we use the gauge conditions in ( 1291 ), we then have 
F grav = -L / N(\a\ 2 + \n\ 2 )d 3 V. 

47r J AH 

The matter flux expression (f2]) of Vaidya solution would be 

F matter : = j^T ab T a £ b Nd 3 V = i- ^ $ 00 Nd 3 V (36) 

where we use 47rT ab r£ & = $g . The total flux of Ashtekar- 
Krishnan then becomes 

Fwa = J [\(t\ 2 + \n\ 2 + <f> 00 ]Nd 3 V. (37) 
Further, the integral of N ^ R can be written as 

[ N {2) Rd 3 V=( dR I {2) Rd 2 V = 87r(i? 2 - Ri) 

J AH J R t J 

where R\ and i?2 are the radii of the horizon at the bound- 
ary cross-sections. For a rotating non-spherical symmetric 
dynamical horizon, we find the relation of the change in the 
horizon area in the dynamical processes can be written as 



dR (R 2 -Rx) 1 



16tt 



N {2) R d 3 V 



AH 



(38) 



2 2 

— -^matter H~ -*grav 

T ab T a t b d 3 V ~\~ ~~ f N(a ab a ab - 2( a ( a )d 3 V. 

AH 1D7T J ah 

Hence, from this equation we can relate the black hole area 
change with energy and angular momentum change. This 
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gives a more general black hole first law in a dynamical space- 
time. If we now define the effective surface gravity 12|] as 



on DH. We summarize our gauge choices and boundary con- 
ditions 



1 

2~R' 



(39) 



then the area of horizon is A = AttR 2 and the differential of 
the area is dA = 8irRdR, therefore 



^dA = -dR. 

8tt 2 



(40) 



For the time evolute vector t a = Nl a - flip a = £° - Qip a , 
the difference of the horizon energy E% can be expressed as 
Jl 



J T ab T a t b d 3 V + yj- / N(* ab <j ab - 2( a C)d 3 V 



1 

16tt 



np ab c^q ab d 3 v. 



AH 



By using ( f3Qb together with the linear combination of 

dR 



T ab T a ^ b d 3 V 



(41) 



(42) 



AH 



-i- f N(a ab a ab - 2U a )d 3 V, 

lOTT J AH 



k = e — e = v = ji — /i = 7 = 7r — a — (8 = 0, 

/9=0,7r=r. (45) 

In order to preserve orthonormal relations, we can choose the 
tetrad as 

£ a = (l,U,X 2 ,X 3 ), n a = (0,-1,0,0), 
m a = (0,0,£ 2 ,£ 3 ). 



in the Bondi coordinate (v, r, x 2 , x 3 ). 

Now we make a coordinate transformation to a new comoving 
coordinate (v, r', x 2 , x 3 ) where r' = r — Ra(v) and Ra(v) 
is radius of a spinning DH. Here 

T = {1 7 U-R A ,X 2 7 X 3 ), n a = (0,-1,0,0), 
m a = (0,0,£ 2 ,£ 3 ), 



where Ra(v) is the rate of changing effective radius of DH. 
From this coordinate, we may see that the dynamical horizon 
is a spacelike or null hypersurface. Here, tangent vector of 
DH is R a = £ a - R A n a =^ where R A > 0. Therefore, it 
implies R a R a < and U, X k = 0(r'). 



we can obtain a generalized black hole first law for dynamical 
horizon 



—dA + nd.J = dE t . 



(43) 



B. The peeling properties and falloff of the Weyl scalars 



Since we use k = v = 0, a ^ 0, A ^ 0, then we have 



m. ASYMPTOTIC EXPANSION FOR A SPINNING 
DYNAMICAL HORIZON 

A. Frame setting and gauge choice 

We choose the incoming null tetrad n a = W a v to be the 
gradient of the null hypersurface v = const. We then have 
9 ' v >a v t a = 0. It gives us the gauge conditions v = /i — ~p = 
7 + 7 = a + /3 — W = 0. Then we further choose n a flag 
plane parallel, it implies 7 = 0. For the setting of outgoing 
null tetrad I, we first choose I to be a geodesic and use null ro- 
tation type III to make e — e = 0. We choose m, m tangent to 
the cross section S, and thus p=p, tt=t 2 . From the boundary 
conditions of a spinning DH ( see eq (HJ), recall 

p=0, tt^O, cr^O. (44) 



2 This also implies uj=0. Kerr solution preferred gauge on horizon is 



(5-4a + 7r)* - (D - 2e)*i=0, (46) 

(A + ^)*o-(^-47f-2/3)1'i=3cr* 2 , (47) 

(D + 2e)cr=* , (48) 

(D + 4e)* 4 - (6 + 4tt + 2a)* 3 = - 3A* 2 , (49) 

(5 + 4/3-r)* 4 -(A + 4^)*3=0, (50) 

(A + 2 A i)A=-1'4. (51) 

in vacuum. Therefore, one can set ^1=^3=0 as peeling prop- 
erties for a spinning DH. This is a similar with perturbation 
method and one may refer to p. 175 and p. 180 in QVD - 

The falloff of the Weyl scalars is algebraically general (this is 
a more general setting than iflill and lfl5ll ) on DH where 

tti = * 3 = 0(r'), * = *2 = *4 = 0(1). (52) 

By considering Vaidya solution as our compared basis for 
matter field part, the falloff of the Ricci spinor components 
are 

$00 = 0(1), $22 = $11 = $02 = $01 = *2i = 0(r')-(53) 
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C. From the radial equations 



(NR10). 



/' 
A 

a 
P 
P 
a 

7T 

r 

e 

u 

k 



X 



*3 



= Mo + (Aio 2 + A A )r' + 0(/ 2 ), 

= Ao + (2/i Ao + *!l)r' + 0(r' 2 ), 

= a Q + [X (W + Pq) + a no}r' + 0(r' 2 ), 

= (3 + [/i 7f + /3 Mo + aoAo]r' + 0(r' 2 ), 

= - 5 7f + TToTfo + croA ]r' + 0(r' 2 ), 

= cr + (-9o7T - ttq + HqVqY + 0(r' 2 ), 

= + [2Mo^o + 2A 7f ]/ + 0(r' 2 ), 

= e + [2a 7f + 2/3 7r + 7f 7r ]r' + 0(r' 2 ), 

= ^° + [Aof°+Moe fe0 ]r' + 0(r' 2 ), 

= 2e r' + 0(r' 2 ), 



2Kr°+^ )/ + 0(r' 2 ). 



*o + (Mo*o - '^°2 + AAV + 0(r' 2 ), 
(-(5 ^ + 37f ^)r' + O(r' 2 ), 
*° + (3^*11 - oto*2V + 0(r' 2 ), 
Mo*2 + (*o-4ft)*V + 0(r /2 ). 



(NR1) 
(NR2) 
(NR3) 
(NR4) 



-i? A (*2 - 9 07T0 + TTOTTO + (T A ) = <J a + $[j 



00- 



- 7T + fTo/io] + 2eo<ro + 

2Ao7fo]_+ 2a TT Q , 
*°] + Sqtto + ttq 



V& 



- Ao(7T 



■A))] 
•A))] 



f cr (o;o 



TTo), 



(TO = i?A[-9o7fo 

Ao = i?A[2/ioAo - 
-2A eo + Mo^o, 
(NR5) d = i?A[aoA < o + 
(NR6) po=RA[a Xo 
(NR7) /t = i?A(Mo + A oA ) + 9 7To 

+(T A - 2/i e + 
(NR8) 5 a = 0, 

(NR9) (5 A - Sofia = Mo^o + A (a - 3/3 ), 
(NR10) Qo^o — QoTfo = 2Im5 7To = — 2Im(Aocro) 
(NR10) Re*!] = -Re[<5 (oo - P )] ~ Re(A cr ) 
+ {a - /3 )(a - Pa)- 

The following equations refer to the equation numbers (a), (b), 
(c), (d) on p. 49 of 01. Here, we relabel (a), (b), (c), (d) as 
(NB1), (NB2), (NB3), (NB4). 



21ml' S], 



(NB1) -5 ^° + (4a - *- )tf° 
= -<5 $8o + ^o$oo. 



Ra[ 



37f * 01 



RA(3fi y° 2 - aa^>l) - Ao*!] + W>$° 



00 • 



(NB2) 

(NB3) -So* 2 - 3tto*^ - fl A (-<5o*2 + fa - 4/3 )*2)> 
(NB4) = 7?a*| - 3A *^ - 4e *i 



$n 
$oi 

$22 
$12 
$02 



= <S>° + 2^ 00 r' + O(r' 2 ), 



RA$i 2 r' 2 + 0(r 
o ji + 0(f /3 

'2 



/3 i 



0(r' 3 ) 
0(r' 2 ) 



D. From the non-radial equations 

The following equations refer to the equation numbers from 
p. 45-p. 47 in H. We re-label (304)-(306) as (NCI), (NC2), 
(NC3). 

(NCI) <5 e = 0, 

(NC2) P = R A [f lQ P + X P}+a P, (54) 
(NC3) P V InP = a -/3 , 



where 

P(v,x k ) :=e 20 



-if , P \7:= So, V:= ^ + »i- 
oar ax 



We relabel (a), (b), (c), (g), (d), (e), (h), (k), (m), (1) as (NR1), 
(NR2), (NR3), (NR4), (NR5), (NR6), (NR7), (NR8), (NR9), 



E. Compatible constant spinor conditions for a rotating 
dynamical horizon 



In this section, we adopt a similar idea of Bramson's asymp- 
totic frame alignment for null infinity |0] and apply it to set 
up spinor frames on the quasi-local horizons. We define the 
spinor frames as 



(55) 



where A^ 

Ai, A as 



Aioa — Aq^a, PA = P\oa — Po L A- We expand 



Ai = X° 1 (v,6,4>) + X{(v,6,4>)r' + 0{r' 2 ), (56) 
Ao = X° (v, 6, <t>) + Aq>, 9, 4>V + 0{r' 2 ). (57) 

Here Ai is type (—1, 0) and Ao is type (1, 0). 

Firstly, we require the frame to be parallely transported along 
the outgoing null normal £ a . 



\im r ,^ o DZ A A =0. 
^rV a (AiOA- A ia) 



0. 



(58) 



Then it gives the condition P Aq =0 on DH. The compatible 
conditions are: 



aoAo* 

9qAi — MoA 



b A° = 0,^ X° a -e X° =0 (59) 

<7 A° = 0, (60) 

= 0,_ (61) 

b A? = -5 A°. (62) 
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IV. SLOW ROTATING BLACK HOLE AND SETTINGS ON 
A TWO SPHERE 



A. Slow rotating Kerr horizon in Bondi coordinate 



Kerr metric in the Eddington-Finkelstein coordinate 

(v,r,6,x) is 



ds 2 



A - a 2 sin 2 i 



S 

2<zsin 2 6>(r 2 + a 2 - A) . a 
— avax + 2a sin 8dxdr 



-dv - 2dvdr 



Y,d8 2 - 



(r 2 + a 2 ) 2 - Aa 2 sin 2 i 



sin 2 9dx ■ (63) 



By changing the coordinate from (i>, 0, x) to (v, 6, x') 



dx = flAdv + dx 



(64) 



where Ha — r i a +a i is the angular velocity on horizon and 
v a is horizon radius of Kerr solution, we can make the term 
g vx dvdx vanished in the 3-D metric. The 3-D metric in the 
new coordinate (v, 8, x') will be 



ds 2 



a 2 sin 2 i 



, , „ a 2 sin 2 _ i , , 
■dv 2 + n^dvdx 



-A 



= • dv 2 - Z A d8 2 - n A 2 



a 2 sin 2 9 , ,2 



■ dv 2 - r 2 A (d8 2 + sin 8 2 d X 



Sa 



(65) 

(66) 
(67) 



and Weyl tensors are: 

K 

P 



7 

7T 

lnV& 2 

#3 
*4 



=0. 



4(t7a - r') 4 
i3\/2a 

J7[(r' - r A ) 2 + a 2 cos 2 + iacosdr') 

2[(r' -r A ) 2 + a 2 cos 2 0] 
[/ _ 1 
2" ~ 47a"' 
ia cos C 



2E A 
a + A 

0,*! =0(r'), 



,7 + 7=0, 



iD cos I 



(r| + a 2 cos 2 



a 2 sin 2 i 



i\/2 sin 8rAf]i 



4S| 



[£>£ 



A + 2ia 2 cos 2 8} 



= 0. 



(71) 
(72) 



a = A = v = 0, 

(ry A - r')(??A - r ') 

- r * + _ r ' ^_I^_J_ (73) 

{r]A - r')(i] A - r') S A r A ' 

iy/2Dif A sin 8 „ iV^DrjA 

r = 



(74) 

(75) 

(76) 

(77) 
(78) 

(79) 

(80) 
(81) 



Remark. In this approximate Kerr tetrad in Bondi coordinate, 
the NP coefficients satisfy 

v = p, — ~p=Tr — a — /3=7 + 7=e — e=0, 

7T = r,p=p,n < 0. (82) 

By examining the approximate Kerr tetrad in Bondi coordi- 
nates, we found it is compatible with our frame setting for the 
asymptotic expansions. 



B. Setting on a two sphere: on horizon cross section 



Here the surface area of slow rotating Kerr is Axerr — 
167rr A . 

Now, we consider the case of slow rotation so that a is small 
and we ignore the a 2 terms. Thus the tetrad components in 
the Bondi coordinate (v, r', 8, (f>) are: 



T = (l,Ur',0,— 
rA 

n a = (0,-1,0,0) 



Dr') (68) 
(69) 

(0,0,1-— ,—(1-—)) (70) 
V2?7a VA sin0 r)A 



To solve the coupling equations from non-radial NP equations 
would be rather complicated and maybe too general to yield 
some interesting physical results. By considering the small- 
tide and slow rotate of DH and consider slow rotate Kerr so- 
lution as a basis from pervious subsection, we use two sphere 
conditions of DH cross section for our later calculation. On a 
sphere with horizon radius Ra(v), one can set 



/'o 



1 

Ra~' 



(83) 



Let P, po on a sphere with radius Ra, then P oc -j^. From 

(NC2), £ fe0 = i?A(M fe0 + Xof°) + vof° which depends 
on the next order nonlinear effect off horizon, we obtain 



where U := ^M. and L> := a(2rA 4 ~ M) ■ The NP coefficients 



Ao 



Ra' 



(84) 
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and 



P = R A noP 



RaP 
Ra ' 



(85) 



Moreover, the effective surface gravity is k = 2cq = jji^' 
and then /i = — 4e (Recall eq. d39l). 

Check the commutation relation [So, Dq]Xo and [So, -Do]co, it 
implies 



R A =0. 



(86) 



By using Komar integral, the quasi-local angular momentum 
on a slow rotating DH is 

J(Ra) = ^-[<f V a <t> b dS ab ]\ A 



— <f> Im( 7T o9o f)dS a (use integration by part) 
~4^/ /Im5 7r d5A 



— jf^ldSA. 



(89) 



This means that the horizon radius will not accelerate (no in- 
flation). The dynamical horizon will increase with a constant 
speed. We note here that if the two sphere condition does not 
hold, then this result is no longer true. 

After applying these conditions, we list the main equations 
that will be used in the later section 



(NR1') 

(NR2') 

(NR3') 
(NR4') 

(NR5') 
(NR6') 
(NR7') 

(NR8') 
(NR9') 

(NB2') 

(NR1')+(NR7') 



RAl-Un + ^2) 



Of) • 



(t = Ra[-^o^o 
+2e a 



*8> 



(-Ra) 2 



a = RacvoHo — oWo, 
$o = RaHo{ttq + /3 ) -+ 
Re* , = 2^ e - 7r 7r 

5 ero = 0, 

— 2oWo = -RaMotto, 
cto9o7to = — 5^AA f o9o7To, 



owo, 
- Re5 7r , 



*° = Ra[^o q 
2ReQ n = |ja 



4J ' 

2fi e . 



a *o + -"o*o 



00 ■ 



From (NB2'), we get Im#P, = sl^ImQoTTo = -3^Im$>°. 
Together with -j^dSA = 2^(1Sa, the angular momentum 



flux for a slow rotating DH is 



J(Ra) 



1 

An 



if 



Ra 
Ra 



f)lm^° 2 dS A 



= -L / Im[(/ - ^f)d ir ]dS A . 



in 



Ra' 



(90) 



We note that from ^ (NR7'), it yields the same result. Here if 

7T ^ and f(v, 9, 0) = G(6, <j))R A (v), then j(R A ) = 0. It 
then returns to the stationary case. If n = 0, i.e., Im\&P> = 0, 
then J and 7 = 0. It then returns to the non-rotating black 
hole. 

For IH, by using integration by parts and Penrose volume I, p. 
281, we have 



For DH, we have 



j. flmdn = f /Im\l/2- 



2W {a R b) rR b 



-W a <j> a 



= -nm a (j) a + C.C. = -ndf + ndf. 



(91) 



(92) 
(93) 



V. ANGULAR MOMENTUM AND ANGULAR 
MOMENTUM FLUX OF A SLOW ROTATING DH 



VI. THE QUASI-LOCAL ENERGY-MOMENTUM AND 
FLUX OF A SLOW ROTATING DH 



Here we use an asymptotically rotating Killing vector 4> a for 
a spinning DH. It coincides with a rotating vector (fi a =ip a 
on a DH and is divergent free. It implies A a </> a := 
SZ S™V b0 (t) a0 = 0. Therefore, 

m a Sr = -mM a . (87) 

Let <fi a = Am a + Bm a , we get A = —B. Therefore, it exists 
a function / such that 

4> a = Sfm a - 5fm a , (88) 

which is type (0, 0). Since / is type (0, 0), therefore 8f = 5/. 



A. Mass and mass flux from Komar integral 



The asymptotic time Killing vector on a DH can be expressed 



as t 



<) 



— dv 



[f a + (U — R A )n a ] | a in corotating coordinate. 



The Komar mass on a DH is then 



Ma 



l - £ W a t h NdS ab = eoNdSt 



4n~f 4Ra~ 



1 



NdS A , 



(94) 



where eo = — Mo/4- This yield the same with eq. (i) from two 
spinor calculation when one chose TV = AqA ,. 
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We then obtain the mass flux on a DH from Komar integral is 
1 



and later we shall see that it agrees with eq. 
spinor method. 



(95) 
from two 



From (ii) We first apply d/dvon (NR7) to get 

*2 + ^otto + Sqttq + 2/3 ir + 2/3 7T /i £o 



Mo = 



2e 



eo 



(98) 



Now, apply time derivative on (ii) and use Sec. IIV Bl vields 



B. Mass and mass flux from two spinor method 



By using the compatible constant spinor conditions for a 
spinning dynamical horizon d60l i, d6TT l and the results of 
the asymptotic expansion, we get the quasi-local energy- 
momentum integral on a slow rotating dynamical horizon 



I(Ra) = -±$H0>%\ydS A 



(i) 

2/3 7ro]AOA°^ A . (ii) 



(99) 

where the total energy momentum flux F tota i is LHS of d99b 
and is equal to matter flux plus gravitational flux F to tai = 
-Pmatter + F gmv . We can write the gravitational flux equal to 
the shear flux plus angular momentum flux. 



F„ 



F a + Fj 



(100) 



In order to calculate flux we need the time related condition 
( 1591 of constant spinor of dynamical horizon in Section ITlI El 
and re-scale it. Then Aq = 0. It's tedious but straightfor- 
ward to calculate the flux expression. It largely depends on 
the non-radial NP equations and the second order NP coeffi- 
cients. By using Sec IIV Bl we substitute them back into the 
energy-momentum flux formula to simplify our expression. 

From (i) Apply time derivative to (i), and then we obtain the 
quasi-local energy momentum flux for dynamical horizon 



(96) 



where it is always positive. Here {iq is the news function of 
DH that always has mass gain. 



From the choice of /i = — we have fi 
^ ( 2 ^R where the two scalar curvature is ^R = -A- (The 

metric of a two sphere with radius Ra is ds 2 = —R\(d0 2 + 
sin 2 8dcj) 2 ).). Integrate the above equation with respect to v 
and use /to = Ra^R/2, we then have (T~ 



d!(R A ) 



1 

8^ 



( 2) R\° XdS A dRA- 



where the shear flux F a is second term of RHS of d99l ) and 
the angular momentum flux Fj is third term of RHS of d99l ). 
The coupling of the shear erg an d ""o can be transform into 
7To terms by using (NR9'). Then integrate the above equation 
with respect to v, we have 



dI(R A ) 



Ra 



}\° X>>dS A dRA. 



(101) 



where dv 



dR A 



Here we note that if one wants to observe 



positive shear flux — ln(i? A croO'o) > 0, it implies that 



°o < 0, 



(102) 



where Ra, Ra > have been considered. So the shear on a 
spinning DH is monotonically decreasing with respect to v. 

Recall that the total flux of Ashtekar-Krishnan d37| i, we com- 
pare our expression with Ashtekar's expression. If we choose 



(97) N = AqA^,, then ( [ToTT i together with d97j gives 



d!(R A ) = 5- 



(2 



RNdS A dR A = ^ / ^{^*oo 
on J R A Ra 



2k 



~Ra~ 



Ra 



}NdS A dR A 



r 



(103) 



where we define \n(R A aoao) := — k for the convenience. 
This is the relation between the change in DH area (Recall 
(l38l ) and total flux with including the matter flux and gravita- 
tional flux. 



Shear flux: In the special case ^ ln(i? A cro< J o) 
k is a constant, we then have 



-Ra^oOo = Ae 



-kv 



-k where 



(104) 
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If k > 0, do \. If k < 0, (To /*■ Therefore, if we want to 
get positive gravitational flux, the shear Co must decrease with 
time v and k > 0. On the contrary, the negative gravitational 
flux implies the shear must grow with time. The negative mass 
loss from shear flux will make the dynamical horizon grow 
with time is physically unreasonable. Therefore, the second 
term of RHS in eq. ( f99b should be positive. This says that the 
shear on a spinning DH will decay to zero when time v goes 
to infinity and the amount of shear flux F a is finite. 



(T -> 0, \v\ -> oo, 
Fu < oo. 



(105) 
(106) 



Hence a slow rotating dynamical horizon will settle down to 
an equilibrium state, i.e., isolated horizon at late time. 

Discussion 

1. If ttq = and the shear does not vanishes (To ^ we 
have 

HRa) = ^- I ^[$oo + 2ka Wo}NdS A . 

47T J R A 

This goes back to the result of the flux of the non- 
rotating dynamical horizon. When k = \, it goes back 
to the result of non-rotating DH in lfl5ll . 

2. If both shear and ttq vanishes, we have 

This result can be compared with dynamical horizon of 
Vaidya solution. 

3. Though we chose the cross section of DH to be two 
sphere, however, it still imply that the shear term can- 
not make into zero. This is because the contribution of 
shear comes from the next order nonlinear effect of the 
equations. 



Laws of black hole dynamics 

LHS of eq. (1103b can be written as 



dI(R A ) 



8tt 2 



(107) 



where A = 4irR A . For a time evolute vector t a = N£ a — 
£lcj) a , the difference of horizon energy dE l can be calculated 
as follow 



16tt7 R a X R a 00 



2k^}N 
Ra 



+ [3A^ tFo - 4-^-Im[(S ./ - ^of)MdV 
Ra -Ka 

and the generalized black hole first law for a slow rotating 
dynamical horizon is 



VII. CONCLUSIONS 



Since ^4 are gauge invariant quantities in a linear pertur- 
bation theory, it allows us to chose a gauge, in which ^3 
vanish on DH. This choice of gauge is crucial for the coupling 
of the NP equations and the consequence physical interpreta- 
tion. In this paper, we use a different peeling property from 
our earlier work lfl5ll ifTill . This leads to a physical picture that 
captures a collapsing slow rotating star and formation of a dy- 
namical horizon that finally settle down to an isolated horizon 
at late time. Further from the peeling property, if the shear 
flux is positive, it excludes the possibility for a slow rotating 
DH to absorb the gravitational radiation from nearby gravita- 
tional sources. The mass and momentum are carried in by the 
incoming gravitational wave and cross into dynamical hori- 
zon. We shall see that though it may exist outgoing wave on 
horizon, however, it will not change the boundary condition or 
make the contribution to the energy flux. A dynamical horizon 
forms inside the star and eat up all the incoming wave when it 
reaches the equilibrium state, i.e., isolated horizon. 

The NP equations are simplified by using two-sphere condi- 
tions for a slow rotating DH with small tide. By using the 
compatibility of the coupling NP equations and the asymp- 
totic constant spinors, the energy flux that cross into a slow 
rotating DH should be positive. The mass gain of a slow ro- 
tating DH can be quantitatively written as matter flux, shear 
flux and angular momentum flux. Further, a result comes out 
that the shear flux must be positive implies the shear must 
monotonically decay with respect to time. This is physically 
reasonable since black hole cannot eat infinite amount of grav- 
itational energy when there is no other gravitational source 
near a slow rotating DH. We further found that the mass and 
mass flux based on Komar integral can yield the same result. 
Therefore, our results are unlikely expression dependent. For 
other quasi-local expressions remain the open question for the 
future study. It would be interesting if one can free the two- 
sphere conditions, then obtain the metric distorted by gravita- 
tional wave. 



-dA + fid J = dE l 



(108) 
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